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Abstract 

Traffic  signals  can  be  synchronized  so  that  a  car,  starting  at 
one  end  of  a  main  artery  and  traveling  at  preassigned  speeds,  can 
go  to  the  other  end  without  stopping  for  a  red  light.   The  portion 
of  a  signal  cycle  for  which  this  is  possible  is  called  the  band- 
width for  that  direction.   Ordinarily  the  bandwidth  in  each  direc- 
tion is  single,  i.e.,  is  not  split  into  two  or  more  intervals  within 
a  cycle.   For  this  case  we  formulate  the  arterial  problem  as  a 
mixed-integer  linear  program:   Given  (1)  an  arbitrary  number  of 
signals,  (2)  the  fraction  of  the  cycle  that  is  red  at  each  signal, 
(3)  upper  and  lower  limits  on  each  signal  period,  (4)  upper  and  lower 
limits  on  speed  each  way  between  adjacent  signals,  (5)  limits  on 
change  in  speed,  and  (6)  a  constant  of  proportionality  between  the 
two  bandwidths,  find  (1)  a  common  signal  period,  (2)  speeds  between 
signals,  and  (3)  the  relative  phasing  of  the  signals  so  as  to  maxi- 
mize the  sum  of  the  bandwidths.   Several  variants  of  the  problem  are 
also  formulated. 

A  branch  and  bound  algorithm  is  developed  for  solving  the  given 
mixed-integer  linear  program  by  solving  a  sequence  of  ordinary  linear 
programs.   A  10  signal  example  is  worked  out. 

The  problem  of  synchronizing  a  network  of  signals  is  also  formu- 
lated.  The  resulting  program  consists  of  the  arterial  programs  for 
the  individual  streets  plus  a  set  of  further  constraints   that  arise 
because  the  streets  connect  together  to  form  closed  loops.   The  objec- 
tive function  used  for  the  network  program  is  a  weighted  sum  of  the 
bandwidths  in  each  direction  on  each  artery.   A  7  signal  example  is 
worked  out. 


1 .   Introduction 

Traffic  signals  can  be  synchronized  so  that  a  car  can  go  from 
one  end  of  a  street  to  the  other  without  stopping,  if  the  driver 
maintains  the  speeds  used  in  setting  the  signals.   The  portion 
of  the  signal  cycle  for  which  this  is  possible  is  called  the  band- 
VJ^dth  for  that  direction.   Synchronizations  for  large  bandwidth  are 
frequently  called  progressions.   Traffic  engineers  have  long  set  up 
progressions  and  they  seem  to  be  quite  helpful  when  traffic  is 
light,  although  possibly  of  little  value  when  traffic  is  heavy. 
The  work  of  Helly  and  Baker     tends  to  suggest  that  this  is  the 
case.   The  value  of  a  progression  can  be  enhanced  by  giving  the 

drivers  speed  instructions,  as  has  been  done  by  von  Stein. 

[3] 
In  an  earlier  paper     j.  Morgan  and  the  author  have  developed 

an  algorithm  for  solving  the  following  two  problems  of  synchronizing 

traffic  signals  for  a  progression  on  an  arterial  street,  in  the 

case  where  bandwidths  in  each  direction  are  unsplit, 

1.  Given  an  arbitrary  number  of  signals,  a  common  signal 
periodj,  the  green  and  red  times  for  each  signal,  and  speci- 
fied travel  times  between  adjacent  signals,  synchronize 
the  signals  to  produce  bandwidths  that  are  equal  in  each 
direction  and  as  large  as  possible. 

2.  Adjust  the  synchronization  to  increase  one  bandwidth  to 
some  specified,  feasible  value  and  maintain  the  other  as 
large  as  is  then  possible. 

Subsequently,  R,  Oliver  pointed  out  to  the  author  that  a  two-signal 
version  of  Problem  1  could  be  set  up  as  a  linear  program.   We  here 
expand  this  idea  into  a  rather  general  formulation  of  the  maximal 
bandwidth  problem.   In  the  general  case  we  still  have  a  linear 


program  but  unfortunately  it  is  of  the  mixed-integer  type.   For 
Problems  1  and  2  the  new  formulation  offers  no  advantages  and  many 
disadvantages.  However,  the  linear  programming  format  opens  up  the 
possibility  of  solving  more  general  problems  involving  the  intro- 
duction of  new  decision  variables  and  new  constraints. 

For  example,  maximal  bandwidth  calculations  as  usually  performed 
have  a  disconcerting  feature.   On  a  long  street  the  critical  signals 
that  constrict  bandwidth  may  turn  out  to  be  very  far  apart.   Then 
a  small  change  in  one  of  the  design  speeds  along  the  street  is 
likely,  upon  re-solving  the  problem,  to  result  in  a  different  syn- 
chronization and  a  different  bandwidth.   Yet,  drivers  do  not  hold 
their  speeds  exactly  constant  and,  indeed,  as  is  well  known,  tend 
to  adjust  their  speeds  to  the  signals.   Thus,  it  makes  a  good  deal 
of  sense  to  let  the  design  speed  between  signals  be  a  variable,  at 
least,  within  certain  limits.   This  can  easily  be  done  in  the  linear 
program. 

Another  variable  that  can  be  introduced  explicitly  is  the 
signal  period.   In  Problem  1,  period  is  a  constant  and,  although 
it  is  not  too  difficult  to  use  our  earlier  methods  to  examine  a 
considerable  number  of  values  in  an  organized  way,  the  continuous 
variation  of  the  linear  program  formulation  seems  preferable. 

Perhaps  the  most  interesting  development  is  that  the  problem 
of  synchronizing  signals  for  a  network  of  streets  can  be  formulated 
as  a  mixed-integer  linear  program.   The  program  for  a  network  con- 
sists of  arterial  programs  for  individual  streets  plus  additional 
constraints  whenever  the  arteries  connect  together  to  form  loops  or 
cycles • 
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2.   The  Arterial  Problem 

2.1   Definitions.   Consider  a  two-way  street  having  n  traffic 

signals.   Directions  on  the  street  will  be  identified  as  outbound 

and  inbound.   The  signals  will  be  denoted  S,^  S„^...,S   with  the 

12      n 

subscript  increasing  in  the  outbound  direction. 

Figure  1  shows  a  space-time  diagram  for  travel  on  a  street. 
Heavy  horizontal  lines  indicate  when  the  signals  are  red.   The  zig- 
zag lines  represent  trajectories  of  the  cars  passing  unimpeded  along 
the  street  in  the  directions  shown.   Changes  in  slope  correspond  to 
changes  in  speed.   The  set  of  possible  unimpeded  trajectories  in  a 
given  direction  forms  a  green  band  whose  horizontal  width  is  the 
bandwidth  for  that  direction.   Although  drawn  but  once^  the  green 
bands  occur  once  per  cycle  in  parallel  bands  across  the  diagram. 

Ordinarily  the  bandwidth  in  each  cycle  is  single_,  i.e.^  is  not 
split  into  two  or  more  intervals  within  a  cycle.   A  split  bandwidth 
can,  however,  occur  and  special  examples  have  been  constructed  in 
which  the  maximal  total  bandwidth  is  made  up  of  two  pieces  in  one  or  both 
directions.     This  possibility  will  be  ignored  in  the  present  formu- 
lation..  The   mathematical  program  to  be  constructed  will  maximize 
the  weighted  sum  of  two  bandwidths,  one  taken  from  ^ach  direction 
without  considering  the  possibility  that  other  pieces  may  exist. 

Certain  signals  form  the  ultimate  limitation  of  bandwidth  and 
will  be  called  critical  signals.   A  signal  S.  is  said  to  be  a 
critical  signal  if  one  side  of  S.'s  red  touches  the  green  band  in 
one  direction  and  the  other  side  touches  the  green  band  in  the 
other  direction.   Thus,  in  Figure  1  signals  S  and  S.  are  critical^ 
but  no  others  are. 
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Figure  1.   Space-time  diagram  showing  outbound  and  inbound  green  bands, 
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2.3  Basic  Maximal  Bandwidth  Formulation.   First  we  set  up 
Problem  1  above  as  a  mixed-integer  linear  problem.   Let 

r.   =   red  time  of  S.  on  street  under  study,  (cycles) 

b  (b)   =   outbound  (inbound)  bandwidth,  (cycles) 
t(h,i)  (t(hji))   =   travel  time  from  S,  to  S.  in  outbound  direc- 
tion (travel  time  from  S.  to  S,  in  inbound 

1     h 

direction) .   These  quantities  are  algebraically 

positive  if  S   follows  S   in  the  outbound 
i         h 

direction,  otherwise  negative,  (cycles) 

0(h,i)  (0(h,i))   =  distance  from  center  of  red  at  S,  to  the 

center  of  a  particular  red  at  S.-   See  Figure 

2.   The  two  reds  are  chosen  so  that  each  is 

immediately  to  the  left  (right)  of  the  same 

outbound  (inbound)  green  band.   0(h,  i)  ('/(h.i)) 

is  positive  if  S  's  center  of  red  lies  to  the 
i 

right  (left)  of  S  's,  (cycles) 
w.  (w.)   =  distance  from  the  right  (left)  side  of  S.'s 

red  to  the  green  band.   See  Figure  2.  (cycles) 
m(h,i)    =  0(h,i)  +  0(h,i)  (2.1) 

Note  that  a  quantity  having  the  dimensions  of  time  can  always  be 
expressed  in  cycles  by  dividing  by  the  period. 
From  Figure  2  can  be  read  the  identities: 


(l/2)r^  +  w^  +  t(h,i)  -  w.  -  (l/2)r.   =  0(h,i)     (2.2a) 
ti    n  1         1 


(l/2)r,  +  w^  +  t(h,i)  -  w.  -  (l/2)r.   =  0(h,i)  .    (2.2b) 
^    '  h    h  1         1 


-4a 


t 

1 


Figure  2.   Geometry  of  green  bands.   Notice  that  0(h^i)  +  0(h,i)  must  equal 
an  integer. 
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The  values  of  0(h, i)  and  0(h,i)  have  the  very  important  constraint 
that  their  sum  must  be  an  integer,  but  otherwise  are  unrestricted. 
Therefore,  the  above  two  constraints  can  be  replaced  by 

(w^+  w^)  -  (w^+  w^)  +  t(h,i)  +  t(h,i)   =  m(h,i)  -  (r  -  r.) 

(2.3) 
m(h, i)   =  integer  . 

From  physical  considerations,  it  follows  that 

t(h,j)   =  t(h,i)  +  t(i,j)  (2.4a) 

t(h,i)   =   -  t(i,h)  (2.5a) 

Therefore, 

0(h,j)   =  0(h,i)  +  0(i,j)  (2.4b) 

0(h,k)   =   -  0(i,h)  (2.5b) 

The  corresponding  inbound  expressions  are  obtained  by  adding  bars 
to  (2.4)  and  (2.5)  above.   Then, 

m(h,j)   =  m(h,i)  +  m(i,j)  (2.4c) 

m(h,i)   =   -  ra(i,h)   .  (2.5c) 

From  Figure  2  we  also  see  that 

w,  +  b  <  1  -  r^  (2.6a) 

w^  +  b  <  1  -  r^  .  (2.6b) 


So  far,  the  expressions  apply  to  any  arbitrarily  indexed 
signals.   No  use  has  been  made  of  the  special  numbering  along  the 
outbound  direction.   In  later  work  on  networks  the  present  gener- 
ality will  be  needed,  but  we  can  simplify  notation  here  by  de- 
fining 

x^  =  x(i,  i-f-1)  (2.7) 

for  X  =   t,  t,  0,  0,  m  . 

Problem  1  will  now  be  represented  by  a  mixed-integer  linear 
program.   Constraints  (2.6)  are  needed  for  each  signal  and  (2.3) 
for  each  adjacent  pair.   (Although  there  are  many  other  types  of 
pairs,  it  suffices  to  consider  adjacent  pairs  since  the  constraint 
(2.3)  for  any  other  pair  can  be  obtained  from  linear  combinations 
of  (2.3)  for  adjacent  pairs.) 

LPl.   Find  b,  b,  w.,  w  ,  m   to 
<^s^^         '   •■   1    i-^   i 

max  (b  +  b)  , 
Subject  to: 


b  =  b 


(LPl.l) 


w.  +  b  <  1  -  r.   ^  (LPl. 2a) 

1  =  1,  .  .  . ,  n 

^i"*"^  ^  ^  -  '^i   ')  (LPl. 2b) 


(W.4.  w.)  -  (w.^^+  w.^p  +  (t.+  t.)   =  m.  -  (r.-  r   )   (lpi.3) 


i+1' 

i=l, . . . , n-1 


ra.   =   integer 
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b,  b,  w,.  w.   >  0 

LPl  has  3n  equations  and  3n-f-  1  unknowns^  not  counting  slack  or 
artificial  variables. 

2,4  Adding  Decision  Variables.   Next  we  let  period  and  speed 
be  variables.   Each  will  be  constrained  by  upper  and  lower  limits. 
In  addition,  changes  in  speed  from  one  street  segment  to  the  next 
will  be  limited.   Finally,  instead  of  requiring  equal  bandwidths 
in  each  direction,  we  let  one  bandwidth  be  a  fixed  proportion  of 
the  other.   (However,  this  is  just  one  of  the  possible  constrain- 
ing relations  that  might  be  used.) 

Let 

k  =  constant  of  proportionality  between  b  and  b. 

T  =  signal  period,  (seconds) 

T, ,  T  =  lower  and  upper  limits  on  period,  i.e., 

T  <  T  <  T   .   (seconds) 

z  =  signal  frequency,  (cycles/second) 

d(h,i)  =  distance  between  S,  and  S..  (meters) 

d.  =  d(i,  i+1) 

V.  (v.)   =  speed  between  S.  and  S.  ,  outbound  (S.  ,  and 
11       '^  1      1+1  1+1 

S   inbound) .   (meters/second) 
i 

e.,f   (e.,f.)  =   lower  and  upper  limits  on  outbound  (inbound) 

speed,  (meters/second) 

1/h  ,l/g.(l/H  ,1/g.)  =  lower  and  upper  limits  on  change  in  outbound 

(inbound)  reciprocal  speed,  i.e. 

1/h.   <   (l/v.^J  -  (1/v.)   <  1/g.  . 
1  —      1+1        1   ~     1 

(meters/second) 
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We  are  constraining  change  in  speed  by  putting  upper  and  lower 
limits  on  change  in  reciprocal  speed.   Although  the  two  are  not 
quite  the  same^  constraining  the  change  in  reciprocal  speed  surely 
satisfies  the  basic  intention,  which  is  to  have  a  means  of  pre- 
venting large  abrupt  speed  changes.   The  reason  for  working  with 
reciprocal  speed  is  that  it  appears  linearly  in  the  constraints  and 

can  be  transformed  into  t    Thus 

i 

t.   =   (d./v.)z  ,      t.  =  (d./v.)z  .  (2.8) 

In  an  expanded  program  t^,t^,  and  z  are  decision  variables,  which, 
once  known,  determine  the  speeds. 

After  some  algebraic  manipulations  the  constrainsts  on  band- 
width, period,  speed,  and  change  in  reciprocal  speed  can  be  added 
to  LPl  to  give  a  more  versatile  mixed-integer  linear  program. 

LP2.   Find  b,  b,  z,  w.,  w.,  t  ,  t"  ,  m   to 
^s^~^~  1   1   i   i   i 

max  (b  +  b) 
subject  to: 


b  =  k  b 


(LP2.1) 


l/T^  <  z  <   1/T^  (LP2.2) 

w.  +  b  <  1  -  r .   ^  (LP2.3a) 


i  =  1,..., 


w.  +  b  <  1  -  r.  J  (LP2.3b) 

(w.+w.)-(w.^^+w.^^)+(t.+r.)  =  m^-(r.-r.^j)\            (LP2.4) 

"i  =  '"''^"  /  i=l,...,n-l 

(d./f.)2  <  t.   <   (d./e^)z  f                               (LP2.5a) 

(d./f.)2  <  t .   <   (d./e.)z  1              (LP2.5b) 


(d./h^z  <  (d./d._^^)t.^^  -  t.  <  (d./g.)2   j  (LP2.6a) 

i=l,...,n-2 

(d./h.)z  <  (d./d._Jt._^  -  t.  <  (d./g.)z  (  (LP2.6b) 

11   —   1   1+1   1+1     1  —   1   1     ' 


b,  b,  w.,  w.   >  0 
'   '   1    1  — 


LP2  involves  (lln-10)  constraints  and  5n  variables,  not 
counting  slack  and  artificial  variables. 

2,5   Symmetric  Problems.   Certain  symmetries,  if  they  exist, 
can  be  used  to  reduce  the  size  of  LP2.   In  particular,  if  the  con- 
straints on  speed  are  the  same  in  each  direction,  and  if  bandwidths 
in  each  direction  are  required  to  be  equal,  the  size  of  the  linear 
program  is  cut  about  407o. 

A  time-saving  technique  may  be  to  solve  for  maximal  equal 
bandwidths  and  then  adjust  the  synchronization  for  other  cases. 
Reference  3  shows  that,  with  some  qualifications  to  cover  special 
conditions,  max  (b  +  b)  is  a  constant  that  can  be  divided  between 
the  directions.   Rules  are  given  there  for  setting  one  bandwidth 
to  any  feasible  value  and  the  other  to  the  largest  value  that  is 
then  possible  for  the  given  speeds  and  period. 

The  linear  program  for  the  symmetric  case  is  based  on  the 
following  two  theorems. 

Theorem  1.   If  the  limits  on  speed  and  change  in  speed  are 
the  same  in  each  direction,  and  if  LP2  has  any  optimal  solutions, 

then  LP2  has  an  optimal  solution  in  which  t.  =  t.  for  each  i. 

1     1 

Proof.   To  have  the  same  limits  on  speed  and  change  in  speed 

in  each  direction  is  to  have  e .  =  e. ,  f .  =  f  ,  e .  =  g  ,  and  h   =  h  , 

1     1    1     1-"  ''i     i'       i     i 
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Let  {t . ,  t,|  be  the  set  of  of  t's  in  an  optimal  solution  for  LP2. 
Define 

t.'  =  t.'  =  (t  +  t.)/2  ,    i  =  l,...,n-l  . 
XI      i    1 

We  claim  that  the  t.  ^t.'   and  t.  <r   t.'   will  yield  a  new  optimal 

11        11        ■'  ^ 

solution  for  LP2.   Since  the  change  of  t's  does  not  affect  the  ob- 
jective function^  the  only  question  is  feasibility.   Clearly, 
LP2.4  is  still  satisfied.   Adding  LP2.5a  to  LP2.5b  and  dividing  by 
2  gives 

(d./f .)z  <  t.'   <   (d,/e.)z  , 
11    —   1   —    11 

so  that  LP2.5a  and  b  are  satisfied.  A  similar  argument  shows  that 
LP2.6a  and  b  are  satisfied.   This  completes  the  proof. 

Theorem  2.   If  the  bandwidths  in  each  direction  are  required 

to  be  equal,  and  if  LP2  has  any  optimal  solutions,  then  LP2  has  an 

optimal  solution  in  which  w.  =  w.  for  each  i. 
1    1 

Proof.   Given  any  optimal  solution  to  LP2  with  b  =  b  ,  replace 

w.  and  w.   bv  w  '  and  w  '  as  defined  by 
X       1    •'   1        X 

w. '   =  w. '   =   (w.  +  w.)/2  . 

L         1  IX 

The  same  type  of  arguments  used  to  Theorem  1  demonstrate  the  feas- 
ibility of  the  new  solution  in  LP2.3a  and  b  and  LP2.4.  Optimality 
is  then  immediate  since  the  objective  function  is  unchanged. 

If  the  conditions  of  both  theorems  are  satisfied  we  can  formu- 
late the  following  simplified  mixed-integer  linear  program; 


LP3.   Find  b,  z,  w. ,  t.,  m.   to 
'      '      1  -^   1    1 


max  b 


-11 


Subject  to; 


(I/T2)   <  z  <   (1/T^)  (LP3.1) 


w.  +  b  <  1  -  r.  ^   i=l,...,n     (LP3.2) 

w.-w    +  t.  =  (l/2)m.  -  (1/2)  (r-r)^  (LP3.3) 

1  a.+l    11  i  1+1  / 


m.  =  integer  (    i=l,...,n-l 

(d./f.)z  <  t.   <  (d./ejz         y  (LP3.4) 

^  1   1    -   1  -    1   i 

(d./h.)z  <   (di/di+i)t._^^  -  t.  <  (d^/8i)zl  (LP3.5) 

J  i=l,...n-2 

b,  w.   >  0 
'   1  — 


LP3  has  (6n-5)  constraints  and  3n  variables  not  counting  slack 
and  artificial  variables. 

2.5  Determining  the  Synchronization.   The  linear  programming 
variables  determine  the  synchronization  of  the  signals.   Let 

e(h^i)   =  relative  phase  (offset)  of  S,  and  S.,  measured 

as  the  time  from  the  center  of  a  red  of  S,  to 

h 

the  next  center  of  red  of  S..  (cycles) 

An  illustration  of  e(h, i)  is  given  in  Figure  2.  We  adopt  the  con- 
vention 0  <  e(h,i)  <  1.  A  set  of  9(1, i),  i=2, ...,n  will  be 
called  a  synchronization  of  the  signals. 

In  order  to  give  simple,  explicit  expressions  for  0  and 
certain  other  quantities,  we  define,  for  arbitrary  real  x: 
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int[x]   =   largest  integer  <  x  (2.9) 

tnan[x]   =   x  -  int[x]  (2.10) 

Thus,  for  example,  int  [5.2]  =  5,   man  [5.2]  =  .1,      int  [-1.3]  =  -2, 
man  [-1.3]  =  .7  . 

From  Figure  2  and  (2.2)  we  see  that 

e(h,i)   =  man[0(h,i)]  (2.11) 

i-1 
0(1,  i)   =  raan[w  -w.  +  E  tj^+(l/2)  (r. -r .)  ]  (2.12) 

k=l  ^ 

In  the  symmetric  case  of  LP3,  (2.2)  shows  that  0(h,i)  =  0(h,i), 
whence 

0(h,i)   =   (l/2)m(h,i)  (2.13) 

i-1 
6(1, i)   =  man[(l/2)   Z   m^]  •  (2.14) 

k=l   ^ 

2.6  Limits  on  the  Integer  Variables.   The  integer  variable 
m(h, i)  is  rather  limited  in  the  values  it  can  take  on  without 
causing  an  infeasible  program.   If  most  of  the  infeasible  values 
can  be  excluded  in  advance,  we  should  be  able  to  save  computation 
time  in  the  algorithm  of  the  next  section. 

We  seek  integers  m' (h^ i)  and  m"(h,i)  such  that 
m'(h,i)   <  m(h,i)   <  m",h,i)  , 
From  (2.3) 

m(h,i)  =  (w,+w,)  -  (w.+w.)  +  t(h,i)  +  t(h,i)  +  (r  -r  )  . 

n   n      1   1  n   1 

Since  0  <  w,   <  1  -  r,   for  any  k  and  since  a  similar  relation 
—   k  —       k 

holds  for  \-i.     , 
k  ' 
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t(h,i)+  t(h,i)-(l-rj^)-(l-r.)  <  m(h,i)  <  t(h,i)+  t(h,i)+  (1-r  )+(l-r.) 

Therefore, 

m"(h,i)   =   int[(l-rj^)  +  (1-r^)  +  t  (h,  i)  +  r(h,  i)  ]       (2.15a) 

m'(h,i)   =  -int[(l-r^)  +  (1-r^)  -  t(h,i)  -  t(h,i)]       (2.15b) 

However,  t  and  t  are  usually  variables  and  can  be  bounded  in  terms 
of  the  constants  of  the  problem.   Thus,  in  LP2 

Following  the  notation  of  LP2,  we  obtain: 

m"^  =  int[(l-r^)  +  (1-r.^^)  +  (d./e.T^)  +  (d./e.T^)]    (2.16a) 

,'^=  -int[(l-r.)  +  (l-r._^^)  -  (d./f.T^)  -  (d./I.T^)].   (2.^6b) 


m' 


3.   A  Branch  and  Bound  Algorithm 

The  algorithms  presently  available  for  solving  the  general  mixed- 
integer  linear  programming  problem  are  usually  considered  only 
partially  satisfactory.   Relatively  small  problems  sometimes  take 
an  inordinate  length  of  time  to  solve.   Consequently,  a  specialized 
algorithm  has  been  developed  for  the  maximal  bandwidth  problem  by 
using  branch  and  bound  methods.   These  methods  have  proven  rather 
successful  in  some  other  combinatorial  problems.  ^^'^^ 

^•^  General  Description.   The  basic  idea  of  the  method  is  to 
break  up  the  set  of  all  feasible  solutions  into  smaller  and  smaller 
subsets  and  to  calculate  for  each  an  upper  (lower)  bound  on  the 
objective  function  of  the  best  solution  therein.   The  bounds  guide 
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the  partitioning  of  the  subsets  and  eventually  identify  a  maximizing 
(minimizing)  solution.   When  a  subset  is  found  that  contains  a  single 
solution  whose  objective  function  is  greater  (less)  than  or  equal  to 
the  upper  (lower)  bounds  for  all  other  subsets^  that  solution  is 
optimal.   The  subsets  of  solutions  are  conveniently  represented  as 
the  nodes  on  a  tree  and  the  process  of  partitioning  as  a  branching 
of  the  tree;  hence  the  name,  "branch  and  bound". 

The  convergence  of  the  process  can  be  assured  by  devising  a 
partitioning  procedure  that^  at  worst^  leads  to  eventual  enumeration 
of  all  solutions.   The  computational  efficiency  of  the  process,  how- 
ever, is  very  dependent  on  the  methods  used  to  perform  the  partitions 
and  calculate  the  bounds.   The  successful  applications  to  date  have 
exploited  special  features  of  the  particular  class  of  problems  at 
hand  to  develop  good  bounds. 

A  branch  and  bound  algorithm  for  LP2  (or  LP3)will  now  be  de- 
scribed.  Consider  an  ^-signal  problem  (2  <  r  <  n)  consisting  of 
the  first  r  signals  of  the  n-signal  problem.   Suppose  that  we  pick 
specific  integer  values  for  the  variables  m  , ...,m  _   in  LP2.   LP2  is 
now  an  ordinary  linear  program  and  can  be  solved  in  a  straight- 
forward manner  to  maximize  the  objective  function.   Now  the  set  of 

n-signal  solutions  that  have  the  specific  values  m^^.'-m    is  a 

1     r-1 

subset  of  the  set  of  all  n-signal  solutions.   Furthermore,  the 
maximal  objective  function  for  the  £-signal  problem  having  these 
m-values  forms  an  upper  bound  on  the  objective  function  for  any 
n-signal  solution  with  these  ra-values,  since  adding  more  signals  can 
only  restrict  the  possibilities  for  maximization.   Thus  we  have  a 
method  for  defining  subsets  of  solutions  (by  specifying  m-values) 
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and  for  placing  bounds  on  their  objective  functions  (by  solving 

ordinary  linear  programs) .   The  subsets  can  be  further  partitioned 

by  specifying  values  for  more  of  the  m's. 

3.2  Flow  Diagram.   The  ideas  of  the  previous  section  will  now 

be  expressed  more  formally.   Let 

X  =  a  subset  of  solutions  of  the  n-signal  problem; 

namely^  the  set  of  all  solutions  that  have  the  m- 

values  of  a  particular  £-signal  problem. 

r(X)   =   the  number  of  signals  in  the  problem  defining  X. 

m,  (X),...,m   fX)   =  the  integer  values  of  m.  for  the  r-signal 
1         r-1  1        — 

problem  defining  X. 
X^  X  ...   =   subsets  of  solutions  obtained  by  partitioning  X 
further. 

The  set  X  is  partitioned  into  subsets  x  ,  X   ...  each  defined  by 
a  different  (r+1) -signal  problem.   Any  of  these  (r+d) -signal  prob- 
lems, say  the  one  for  X.,  will  have  the  same  m  's  as  X  up  to 
i  =  r-1: 

m. (X.)   =  m. (X)         i=l, .. .r-1  , 

and.  in  addition,  a  value  for  m  .   Let 

3  r 

m  '   ID  "  =  smallest  and  largest  value  of  m  to  be  used  in 
r    r  "  r 

partitioning  X. 
a(X)   =  an  upper  bound  on  the  objective  function  for  the 

solutions  in  X^  as  obtained  by  solving  an  ordinary 
linear  program  for  the  r-signal  problem  with 
specified  m-values. 
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p  =  best  n-signal  objective  function  found  so  far  in 
the  calculation. 

A  flow  diagram  of  the  algorithm  is  shown  in  Figure  3.   We 
shall  trace  through  the  boxes.   As  this  is  done,  it  may  be  helpful 
to  refer  to  Figure  h,    which  shows  the  tree  of  a  worked-out  10 
signal  problem.   In  the  tree  representation,  subsets  of  solutions 
are  represented  by  nodes  and  so  we  shall  use  X  to  refer  to  both 
a  subset  and  its  node. 

Box  1  initializes  three  quantities:   X  to  the  set  of  all 
solutions,  r   to   1,  and  p  to  zero.   Box  2  initializes  two  quanti- 
ties to  be  used  in  the  construction  of  the  nodes  branching  out  of 
X:   The  index  j  for  counting  the  nodes  is  set  to  1  and  m  ,  whose 

values  will  define  the  nodes,  is  set  to  its  lower  limit  m  '. 

r 

Box  3  sets  up  a  node  X  on  a  branch  out  of  X.   The  node  may 
J 

be  identified  by  m  (X.),  which  is  set  to  the  current  value  of  m  . 

r  J  '  r 

(In  Figure  4  the  value  of  m  (X  )  is  written  inside  the  circle  of 

r   j 

each  node,  x  . •   To  define  the  node  completely,  we  also  need 

m(X. ),..., m  ^(X.),  but  these  can  be  found  by  tracing  back 

through  the  tree  to  "all  solutions"  and  reading  off  the  m.  in  the 

1 

nodes  encountered.   Next  LP2  is  solved  for  the  (r+1) -signal  prob- 
lem defined  by  the  current  m-values,  m, ,...,m  .   The  maximal 

—       '   1      r 

value  of  the  objective  function  becomes  the  upper  bound  a(X.)  and 
is  used  to  label  the  node. 

Box  4  tests  whether  an  n-signal  problem  is  currently  being 
solved.  If  so,  box  5  tests  to  see  whether  the  resulting  objec- 
tive function  is  better  than  the  previous  best.  If  this  is  the 
case,  box  6  replaces  the  previous  best  with  the  new  solution. 


-16a 


-> 


^— > 


i 

t 

S.t:- 

0    ^^ 

[ 

— ; 


v 


LT&p.3lir^- 


V 


<r 


vv-C> 


;"  vH  -I7  I  w^=  wv':?l 


Zl- 


v^^ 


A 


^    }     '^t^-^tl^     \j,-^\     ■t><-y%-\Vv-<Kj.     V-C^X 


i     ^ 


';^\i»^- 


'Aw^ 


V  <  v\     -tor  0.    wc-fiJ- 


-i 


C^  o^   .         SiJ    Z    tb  ^.  ti^ 


4^- 


I 


',10  1     o<  IK)  ■> 

fc, «  -      ^- 


'•n 


Figure  3.   Flow  diagram  of  a  branch  and 
bound  algorithm  for  LP2. 
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Box  7  tests  to  see  whether  all  the  planned  m  values  have 
r 

been  used  to  create  new  nodes.   If  not,  box  8  increases  i  and  m 
by  one  and  returns  to  box  3  to  set  up  a  new  node. 

Box  9  searches  for  a  new  node  from  which  to  branch.   The 
candidates  are  the  terminal  nodes  of  the  tree,  i.e.,  those  with- 
out branches.   There  is  no  need  to  consider  nodes  with  r  =  n. 
Of  the  eligible  nodes,  the  one  with  the  largest  a  is  selected  and 
becomes  the  new  X.   The  current  value  of  r  is  reset  to  r(X)  and 
the  current  values  of  m.  to  m. (X)  for  i=l,,.,,r-l. 

Box  10  tests  to  see  whether  the  largest  upper  bound,  a,  is 
greater  than  the  best  n-signal  objective  function  so  far.   If  so, 
there  is  some  possibility  of  finding  a  larger  n-signal  objective 
function  and  the  calculation  returns  to  box  2  to  branch  from  the 
new  X.   Otherwise,  the  calculation  is  finished,  and  the  current 
best  solution  is  optimal  for  the  original  mixed-integer  n-signal 
problem. 

It  may  be  helpful  to  trace  through  the  first  few  nodes  in 
Figure  4.   First  the  node  "all  solutions"  is  laid  out.   Then, 
since  m  '  =  0  and  m  "  =  1,  the  nodes  m  =  0  with  a  =  .440  and 


m^  =  1  with  a  =  .292  are  laid  out.   The  largest  a  of  these 


two 


is  .440  and  so  the  m  =  0  node  is  for  branching.   Out  of  it  comes 
the  nodes:  m  =  0,  a  =  .  282  and  m  =  1,  a  =  .381.   The  next  step 
would  be  to  branch  from  the  m„  =  1,  a  =  . 381  node.   Eventually 
the  tree  is  developed  as  shown.   Notice  that  the  10-signal  node 
labeled  "optimum"  has  an  objective  function  of  .282,  which  is 
greater  than  or  equal  to  the  upper  bound  for  any  of  the  other 
terminal  nodes. 
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Figure  4.   Tree  for  10  sigi.al 
arterial  problem.   Labels  on 
nodes  are  the  calculated  upper 
bounds  on  bandwidth  (in  cycles) 
for  the  solutions  represented 
by  the  node.   The  label  X  in- 
dicates no  feasible  solution. 
The  numbers  inside  the  nodes 
are  values  for  the  m.  shown  at 
the  right. 
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3.3  Example.   The  10  signal  problem  being  used  as  a  numeri- 
cal example  is  taken  from  Reference  3  and  represents  a  stretch  of 
Euclid  Avenue  in  Cleveland,   The  distances  and  green  splits  are 
taken  from  actual  usage.   The  limits  on  period,  speeds,  and  speed 
changes  have  been  made  up  fairly  arbitrarily. 

The  constants  of  the  problem  are  briefly  as  follows.   The  sig- 
nals, starting  with  S,  and  workipg  outbound,  are  located  at  0,  168, 
381,  716,  929,  1173,  1371,  1493,  1706,  and  1843  meters.   The  cor- 
responding red  times  are  .47,  .40,  .40,  .47,  .48,  .42,  .40,  .40,  .40, 
and  .42  cycles.   Limits  on  period  are  T  =  55  seconds  and  T  =  75 
seconds.   Lower  and  upper  limits  on  speed  are  e.  =  e.  =  13.4  meters/ 
second  (30  mph,  48.3  km/hr)  and  f.  =  f.  =  17.9  meters/second  (40  mph, 
64.4  km/hr)  for  each  i.   Limits  on  change  in  reciprocal  speed  are 
1/g.  =  1/g.  =  "  1/h.  =  r   1/h.  =  .0121  (meters /second)   for  each  i. 
This  corresponds  to  a  maximum  possible  change  in  speed  of  about  ±2.2 
meters/second  (4.9  mph,  7.9  km/hr)  at  the  lower  limit  of  speed  and 
±3.9  meters/second  (8,7  mph,  14.0  km/hr)  at  the  upper  limit  of  speed. 

The  problem  can  be  solved  by  LP2  or,  since  we  have  chosen  to 
solve  for  equal  bandwidths,  by  LP3.   The  tree  developed  by  the  branch 
and  bound  algorithm  appears  in  Figure  4.   A  total  of  52  linear  pro- 
grams are  required.   The  maximal  bandwidth  for  each  direction  is  .282 
cycles.   The  space-time  diagram  for  the  street  is  shown  in  Figure  5. 
The  optimal  period  is  75  seconds,  making  the  bandwidth  21.2  seconds. 
Starting  between  S,  and  S^  and  working  up,  the  speeds  are:   17,9, 
17.9,  17.1,  14.2,  13.4,  14,9,  13.4,  15.6,  17.9  meters/second.   The 
phases,  6(1,  i),  in  order  of  increasing  index  starting  with  i=l, 
are:   0,  0,  0,  1/2,  1/2,  0,  0,  0,  1/2,  1/2  cycles.   The  critical 
signals  are  S  ,  S.,  S  ,  S,,  and  S.  . 
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The  results  can  be  compared  to  those  obtained  by  exploring 
various  constant  speeds  and  periods  with  the  methods  of  Reference 
3.   The  best  result  found  was  15.2  meters/second  at  65  seconds 
for  a  bandwidth  of  .235  cycles.   This  is  essentially  the  syn- 
chronization displayed  in  Reference  3.   The  additional  flexibility 
of  variable  speed  has  permitted  at  207.,  increase  in  bandwidth. 

3.4  Discussion.   For  the  numerical  example  reported  here^  the 
steps  of  the  branch  and  bound  algorithmwere  carried  out  manually 
except  for  the  linear  programs.   In  some  linear  programming  codes^ 
it  is  possible  to  suppress  constraints.   Then  the  entire  n-signal 
problem  can  be  loaded  and  only  those  constraints  corresponding  to 
the  current  £-signal  problem  used.   Changes  in  the  m.  are  ef- 
fected by  changes  in  the  constants  vector.   Frequently^  time  can 
be  saved  by  starting  a  new  problem  from  the  basis  of  a  similar 
problem  just  completed. 

The  computational  limits  of  the  algorithm  are  relatively  un- 
explored.  If  experience  from  a  rather  different  application 
is  a  guide^  the  number  of  linear  programs  required  is  likely  to 
increase  exponentially  with  n.   If  so^  there  may  well  be  a  fairly 
sharp  upper  limit  to  the  size  of  the  problem  that  is  computation- 
ally feasible.   We  feel^  however,  that  the  example  here  demon- 
strates that  problems  of  a  size  of  practical  interest  can  be  solved. 

A  solution  that  is  likely  to  be  good  but  may  not  be  optimal 
can  be  obtained  by  solving  roughly  2n  linear  programs.  The  pro- 
cedure is  simply  to  branch  next  from  the  node  with  the  largest  a 
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among  those  just  created  from  the  current  X.   Frequently  there 
are  only  two  such  nodes.   The  scanning  of  the  other  terminal 
nodes  is  omitted.   The  value  of  r  increases  raonotonically  to  n. 
In  the  example  of  Figure  4,  the  solution  obtained  in  this  manner 
has  a  bandwidth  of  .-278. 

4.   Network  Problems 

The  mixed-integer  formulation  can  be  extended  to  networks. 
The  network  program  consists  of  arterial  programs  for  the  in- 
dividual streets  plus  "cycle  constraints"  that  connect  the 
streets  together.   An  objective  function  can  be  formed  from  the 
bandwidths  of  the  arteries.   A  new  decision  variable,  the  red- 
green  split,  can  usefully  be  introduced  at  certain  signals. 

Some  of  the  decision  variables  of  LP2  will  be  left  out  of 
the  network  formulation  to  be  developed  in  this  section.   The 
reason  is  that  we  can  expect  network  problems  to  be  larger  and 
computationally  more  difficult  than  individual  arterial  problems. 
Therefore  it  seems  desirable  to  show  ways  in  which  the  program 
might  be  simplified.   We  shall  keep  speed  fixed  over  the  length 
of  an  artery,   but  variable  both  for  the  artery  as  a  whole  and 
from  artery  to  artery. 

4.1  Arterial  Constraints.   In  a  network  problem  the  term 
artery  will  be  used  to  refer  to  any  street  in  which  a  progression 
is  desired,  i.e.,  on  which  we  define  a  bandwidth  and  bring  it 
into  the  mathematical  program.   As  earlier,  the  signals  of  the 
network  are  presumed  to  be  designated  S^,  S  ,...,S  but  we  can 
no  longer  suppose  that  signals  with  adjacent  index  values  are 
adjacent.   Consider  then  an  artery  and  suppose  that  the  signals 
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at  its  extreme  ends  are  S.  and  S..   We  designate  one  direction, 

^      J 

say  S.  to  S.;  and  outbound  and  thereby  define  artery  ij. 

Quantities  relating  to  artery  ij  will  be  identified  by  a 
superscript  ij.   Thus,  in  a  straightforward  way,  the  appropriate 
definitions  of  Section  2  generalize  to  b  ,  r   ,  and  w    .   In 
addition,  let 

V    =   speed  in  outbound  direction  on  artery  ij, 

(meters/second) 

ij   ij 
e  ,  f    =   lower  and  upper  limits  on  outbound  speed  on 

artery  ij.   (meters/second) 

u    =   1/v  T  =  reciprocal  speed  outbound  on  artery 

-1 
ij.  (meters/cycle) 

Analogous  notation  is  used  as  appropriate  for  inbound  variables. 
When  two  signals  are  required  to  define  a  quantity,  the  arterial 
designation  is  redundant  and  will  be  omitted.   For  example,  the 
quantities  0(h,i),  0(h,i),  m(h, i)  and  d(h, i)  will  be  used  as  de- 
fined in  Section  2. 

The  constraints  for  artery  ij  are  as  follows.   For  the  artery 
as  a  whole: 


(l/f^^)z  <  u^^  <  (l/e'"^)z   ,  (LP4.1a) 

(l/f'"-^)z  <  u^^  <   (l/e''-^)z   ;  (LP4,lb) 


for  each  signal  S  on  the  artery: 

Wj^^^  +  b^^  <  1  -  t^^^  ,  (LP4.2a) 

w^^^     +  b^-^  <  1  -  r^^^  ;  aP4.2b) 
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for  each  pair  of  adjacent  siganls  on  the  artery,  say,  S   followed 


by  S „  outbound: 


(Wj^^^+  Wj^""^)  -  (w/J+  w/J)  +  d(k,Q)  (u'-J+  u'-J) 

(LP4.3) 

m(k,Jl)      =   integer  (LP4.4) 


4,2  Objective  Function.   In  a  progression  the  bandwidth  de- 
termines the  maximum  length  (in  time)  of  a  platoon  that  can  pass 
unimpeded  down  the  street  when  the  prescribed  speeds  are  main- 
tained.  Therefore,  when  a  progression  system  is  appropriate,  the 
streets  (or  directions  on  streets)  should  presumably  be  allocated 
bandwidth  in  some  manner  according  to  the  flows  on  them.   Two 
devices  are  available  for  allocating  bandwidth  among  streets  (and 
directions):   the  objective  function  and  the  constraints. 

Ideally  we  might  constrain  each  bandwidth  to  be  greater  than 
some  specified  size  needed  to  pass  a  known  flow.   Unfortunately, 
there  is  no  guarantee  that  the  result  will  be  feasible.   As  an 
alternative,  we  can  choose  some  important  bandwidth,  call  it  b  , 
and  maximize  this,  but  require  that  each  other  bandwidth  be 
greater  than  or  equal  to  a  specified  fraction  of  b  .   The  frac- 
tions would  be  chosen  according  to  relative  flows. 

Such  constraints  are  not  sufficient  to  get  the  most  out  of 
the  system.   Sometimes  the  bandwidth  on  an  artery  will  not  be  in 
conflict  with  bandwidth  on  other  arteries.   In  order  to  guarantee 
that  such  a  bandwidth  is  maximized,  every  bandwidth  must  be  in- 
cluded in  the  objective  function.   The  less  important  bandwidths 
can  be  weighted  by  small  coefficients  and  b   by  a  large  coefficient, 
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o   1     q 

To  formalize  these  ideas^  let  b  ,  b  ^...,b  be  the  band- 
widths  of  the  problem  (reindexed  here  for  convenience)  and  let 

a   =  weight  assigned  to  the  i   bandwidth;  i=0^....q, 

i  o  ,     ,    th 

k   =  fraction  of  b   to  be  guaranteed  to  the  i   band- 
width; i=l, . . . ,  q. 

We  take  as  the  objective  of  the  network  program 


o  o  Q  q 

max  a  b  +  . . ,  +  a^b^  ,  (LP4.5) 


as  constrained  by 


b^  >  k^  b°  ,  i=l, .. .,q  (LP4.6) 

1 

Some  care  must  be  exercised  in  the  choice  of  the  k  so  as  not  to 

o 
let  an  unimportant  artery  limit  b  . 


4.3  Cycle  Constraints.   A  new  constraint  with  a  new  integer 
variable  must  be  introduced  whenever  several  arteries  intersect 
to  form  a  closed  loop  or,  as  we  shall  call  it,  a  cycle.   An  as 
example  of  a  cycle,  S.^  S  ,  S  and  S  appear  at  the  arterial  in- 
tersections of  a  cycle  in  Figure  6.   The  basic  reason  for  the  new 
constraints  lies  in  a  physical  requirement  of  simultaneity. 
Suppose  we  set  a  master  clock  to  zero  at  a  center  of  red  of  S  . 
If  we  proceed  down  artery  13  to  S  ,  the  time  of  the  center  of  red 
there  is  fixed  by  the  progression  on  13.   As  a  result,  the  time 
of  the  center  of  red  for  S„  along  the  direction  of  artery  35  is 
also  fixed  on  the  master  clock.   The  progression  on  35  then 
fixes  the  time  of  center  of  red  at  S^.   Continuing  around  the 
cycle,  we  eventually  return  to  determine  the  time  of  a  center  of 
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red  at  S  .   This  time  must  be  an  integral  number  of  periods  after 
time  zero.   The  algebraic  statement  of  this  requirement  consti- 
tues  the  new  constraint. 

The  terminology  of  graph  theory  is  helpful  in  describing 
networks.   Signals  appear  at  nodes.   Lengths  of  street  between 
adjacent  signals  are  called  arcs .   An  arc  has  a  direction^  which 
will  be  the  direction  chosen  as  outbound  on  the  street  segment. 
A  sequence  of  arcs^  irrespective  of  arc  directions^  that  closes 
on  itself  will  be  a  cycle.   The  node  for  S  will  be  denoted  simply 
by  k.   The  arc  from  i  to  j  will  be  denoted  (i,j).   An  arc's  direc- 
tion is  indicated  in  diagrams  by  an  arrow. 

The  formulation  of  the  constraint  for  a  single  cycle  pro- 
ceeds as  follows:   Suppose  we  trace  around  a  cycle  and  find 

pieces  of  arteries  i,j,,  i«Jo;---i  ^   i"  that  order.   Suppose  that 

1  1   2  -^     p  p 

the  signals  at  which  arteries  intersect  are  represented  by  the 

nodes,  k, ,  k  .....  k  with  k,  at  the  intersection  of  arteries 
'12    '  p       1 

i  j   and  i,j,;  k-  at  the  intersection  of  i.j..  and  i^j„^  etc.   The 
cycle  will  then  be  denoted  C(k. ,...^k  ), 

Let  t  =  0  be  the  time  (measured  in  cycles)  of  a  center  of  red 
at  k  along  artery  i.j^.   Then  t  =  0(k  ,  k  )  is  the  time  of  a 
center  of  red  at  k„  along  artery  i,j,-   The  same  instant  is  also 
the  center  of  gifeen  at  k  along  artery  i„j-.   Therefore  a  center 


2^2' 


of  red  at  k  along  i  j   occurs  at 


0(k^,  k^)  +  1/2.* 


*  The  use  of  1/2  assumes  that  red  time  on  one  street  of  the  inter- 
section coincides  with  green  time  on  the  cross  street  and  vice 
versa.   More  complicated  arrangements  are  possible  and  might  re- 
quire replacing  the  1/2  by  some  other  fraction. 
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Starting  from  this  time  we  can  add  0(k„^k  )  to  find  a  center  of 
red  at  k  along  artery  i  j   .   Adding  1/2  we  turn  the  corner  and 
get  on  artery  ioJ-,  •   Proceeding  around  the  cycle,  we  eventually 
determine  the  time  at  which  a  center  of  red  occurs  at  k  along 
artery  i  j   .   Let 

c(k  , ...,k  )  =  integer  variable  for  the  cycle  C(k  , ...,k  ). 

Then  the  cycle  constraint  is: 


c(k,,...,k  )  =  0(k  k  )+0(k  k  )+...+0(k    k  )+0(k  ,kJ  +  (p/2).   (LP4.7) 
i       p        12      Li  p-1   p      P   1 

Appropriate  expressions  for  0(i,j)  in  terms  of  other  variables  of 
the  program  are  given  by  (2.2)  or  (2.13)^  possibly  with  some 
help  from  (2.4)  and  (2.5) . 

Multiple  cycles  introduce  multiple  constraints.   The  number 
of  possible  cycles  is  greater  than  the  number  of  required  con- 
straints.  The  minimum  number  of  constraints  can  be  developed  by 
tracing  over  the  network  as  follows:   Pick  a  starting  node.   Trace 
over  the  arcs  (irrespective  of  their  direction)  that  make  up  each 
artery  passing  through  the  starting  node.   Then,  for  each  node  reached 
by  the  tracing  so  far,  trace  over  the  arcs  on  the  arteries  that  pass 
through  the  node.   Now,  whenever  the  tracing  procedure  closes  a 
cycle,  form  a  constraint  from  the  cycle.   Continue  the  process 
until  no  further  arcs  can  be  traced.   If  the  whole  network  has 
been  traced,  the  job  is  finished.   Otherwise,  the  network  breaks 
down  into  two  (or  more)  disconnected  sub-networks,  one  of  which 
has  just  been  isolated.   Pick  another  starting  node  and  continue. 

The  tracing  process  works  because  it  identifies  which  signals 
bear  rigid  time  relationships  to  the  starting  signal  by  virtue  of 
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being  connected  to  it  through  a  sequence  of  arterial  progressions. 
Whenever  a  cycle  is  closed  the  rigid  relationships  must  be  made 
compatible  by  a  constraint.. 

If  the  network  has  a  signal  at  every  intersection  of  arteries 
(i.e.  there  are  no  overpasses),  the  graph  of  the  system  is  planar. 
Then  the  number  of  cycle  constraints  equals  the  number  of  distinct 
areas  enclosed  by  arteries.   Figure  6  shows  a  planar  graph  re- 
quiring two  cycle  constraints.   Cycle  constraints  are  not  neces- 
sarily unique.   If  Figure  6,    for  example,  we  could  use  the  pair 
from  C(l,3,4,7)  and  C(7,4,5,6)  or  equally  well  the  pair  from 
C(l,3,4,7)  and  C(l,3,5,6).   Algebraically,  either  of  these  pairs 
can  be  obtained  from  the  other. 

4.4  Red  Time  Constraints.   A  signal  at  the  intersection  of 
two  arteries  may  be  a  critical  signal  for  one  but  not  the  other. 
Then  a  shift  of  green  time  from  a  street  to  the  cross  street  may 
increase  bandwidth  on  one  artery  without  reducing  it  on  the  other. 
Thus,  the  red  time  at  arterial  intersections  may  sometimes  be  a 
useful  decision  variable.   The  variable  will  presumably  be  con- 
strained.  The  fraction  of  time  red  on  each  street  may  be  limited 
because  of  flow  requirements.   The  absolute  time  red  may  be  limited 
on  the  low  side  to  permit  pedestrian  crossing,  and  on  the  high  side 
to  placate  drivers. 

For  signals  at  which  red  time  is  a  decision  variable,  let 

E   ,  F,     =   lower  and  upper  limits  on  r    .  (cycles) 


G,   ,  H,     =   lower  and  upper  limits  on  r   T.  (seconds) 
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Then  the  variable  r    is  constrained  by 

ij       ij       ij 
\    <  \    <  Fj^  (LP4.8a) 


ij        ij       ii 
\^     <     ^^         <     \     z  (LP4.8b) 


4.5   Steps  in  Formulating  a  Network  Problem.   The  general 
mixed-integer  linear  program  for  the  network  problem  is  notation- 
ally  soraevhat  cumbersome  and  rather  than  write  the  program  out  we 
list  the  steps  that  generate  it. 

1.  Determine  which  streets  of  the  network  are  to  be 
arteries, 

2.  Set  up  the  objective  function  and  related  con- 
straintSj  (LP4.5)  and  (LP4.6) 

3.  Set  up  the  period  constraints; 

(L/T2)   <  z  <   (1/Tj^)  (LP4.9) 

4.  For  each  artery,  set  up  the  constraints, 
(LP4.1)  -  (LP4.4),  or,  if  preferred,  the  constraints  of  LP2  or 
LP3. 

5.  Set  up  cycle  constraints  using  (LP4.7)  and  the 
procedure  of  Section  4.3. 

6.  Add  red  time  constraints  (LP4.8)  for  any  red  times 
being  made  decision  variables. 

Network  problems  can  be  solved  by  the  branch  and  bound  methods  of 
Section  3  with  a  few  modifications.   As  there,  the  integer 
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variables  are  introduced  one  at  a  time.  At  any  stage^  specific 
values  are  given  to  some  of  the  integer  variables  and  serve  to 
define  a  subset  of  solutions  to  the  network  problem.  An  upper 
bound  on  the  objective  function  for  solutions  in  the  subset  is 
found  by  solving  an  ordinary  linear  program  using  the  specific 
values  of  the  integer  variables. 

However,  in  the  network  problem,  the  full  n-signal  objective 
function  must  be  used  even  for  an  r-signal  problem.   Otherwise  the 
upper  bound  obtained  is  inappropriate.   As  a  result,  however, 
every  bandwidth  that  appears  in  the  objective  function  must  have 
a  constraining  equation  (LP4.2)  in  each  linear  program  solved,  for 
otherwise  the  objective  function  will  have  an  unbounded  maximum. 
Therefore,  the  "£-signal"  problem  will  usually  include  constraints 
from  more  than  r  signals.   It  would  in  fact  be  permissible  to  use 
all  the  constraints  of  the  original  problem  (except  those  involving 
integer  variables  not  yet  assigned  specific  values)  in  every 
ordinary  linear  program.   With  the  exception  of  the  unboundedness 
problem  just  noted,  however,  the  constraints  associated  with  a 
signal  would  not  be  expected  to  affect  the  objective  function  very 
much  until  the  constraint  (LP4.3)  containing  the  integer  variable 
is  added.   Thus,  adding  the  constraints  signal  by  signal  holds 
down  the  size  of  the  linear  program  and  may  be  expected  to  save 
computation  time. 

The  order  of  adding  signals,  however,  will  affect  the  size  of 
the  tree.   By  and  large,  troublesome  constraints  should  be  intro- 
duced early.   Examples  would  be  constraints  for  signals  close 
together  and  constraints  for  cycles  that  have  small  perimeters. 
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4.6  Example.   Figure  6  displays  a  7  signal,  5  artery  prob- 
lem which  has  been  solved  in  the  symmetric  case.   The  arteries 
are  13,  35,  47,  56,  and  16.   The  objective  function  used  is 

13    ^,,,35   ,16   ,56   ,47, 
b   +  .01(b   +  b   +  b   +  b  ), 

along  with  the  constraints 

,35    „13   ,16    ^,13   ,56    ^,13   ,47    ^,13 
b   >  .5b  ,   b   >  .5b   ,   b   >  .5b   ,   b   >  .5b 

Limits  on  period  are  50  and  100  seconds.   Speed  on  each  artery  is 

limited  by  14  and  16  meters/second  (31.3  and  35.8  mph) .   All  red 

13 
times  are  taken  to  be  .5  cycles  except  r   ,  which  is  made  a  de- 

13 
cision  variable  subject  to:    -^  £  ^7   S  -^     cycles  and 

13 
25  <  r   T  <  50  seconds.   Distances  between  signals  are  shown  in 

Figure  6. 

The  complete  program  is  given  in  the  appendix  along  with  the 
solution  tree.   There  are  44  constraints  and  10  integer  variables. 
However,  in  the  branch  and  bound  solution  we  deal  with  42  con- 
straints and  8  integer  variables.   The  reason  is  that  each  of  the 
two  cycle  constraints  is  introduced  by  limiting  the  values  that 
some  regular  integer  variable  is  permitted  to  range  over. 

The  network  example  is  not  hard  to  solve.   15  linear  programs 

13 
were  required.   In  the  optimal  solution  the  bandwidths  are  b   =  .35, 

b   =  .286,  b   =  .5,  b   =  .5,  and  b   =  .286  cycles.  The  period 

is  62.5  seconds.   The  speeds  are  v   =  16,  v   =  14,  v  =  16, 

V   =  16,  and  v   =14  meters/second.   The  value  of  r  =  .5 

cycles . 
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Figure  6.   A  network  of  7  signals.   In  this  network,  five  arteries:   13, 
35,  56,  16,  and  47,  have  been  selected  for  consideration.  Two  cycle  con- 
straints are  required.   Arrows  indicate  direction  chosen  as  outbound. 
Numbers  on  street  segments  are  distances  in  meters.  The  red-green  split 
has  been  made  a  decision  variable  at  7. 
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5.   Summary 

A  variety  of  maximal  bandwidth  problems  have  been  formulated 
as  mixed-integer  programs.   These  include  single  arterial  problems, 
more  complicated  arterial  problems  using  speed  and  period  as  de- 
cision variables,  and  problems  involving  street  networks.   Some  of 
these  cases,  including  that  of  the  street  network,  do  not  appear  to 
have  been  put  in  any  kind  of  optimizing  format  before. 

Branch  and  bound  methods  are  given  for  solving  the  programs 
generated  and  examples  are  worked  out.   Although  it  is  not  known 
at  this  point  how  large  a  program  can  reasonably  be  solved  by  these 
methods,  it  seems  likely  that  problems  large  enough  to  be  of 
practical  interest  can  be  solved. 
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Appendix  1 

The  raixed-int.eger  linear  program  for  the  network  of  Figure  6 
is  given  below  for  the  symmetric  case. 

Find  b  ,  z,  u  •^,  w^  ,    m(k,J?),  c  (i,j,k,^),  and  r^^   to 


ul3  „wu35        ,16        ,56        ^47, 

max  b        +.01(b        +b        +b        +b) 


Subject  to: 

,^   .  .  .  ,.        .      ,35        13 

bandwidth  constraints:   b    >   . 5b  (A.l) 

b    >   .5b  (A. 2) 

^56      ^,  13 

b    >   .5b  (A. 3) 

,  47        13 

b    >   .5b  (A. 4) 


period  constraints:      .01  <  z  <   .02  (A. 5, 6) 

13 
artery  13:  .0625z  <  u    <   .0714z  (A. 7, 8) 

13    13 
w    +  b    <   .5  (A. 9) 

13    13 
W2   +  b    <   .6  (A. 10) 

"l^"^"  ^2^"^  "•"  200u^^  =   .5m(l,2)  -  .05   (A. 11) 

13    13 
w    +  b    <   .5  (A. 12) 

W2^^-  w^^^  +  300u^"^  =   .5m(2,3)  +.  05    (A. 13) 
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artery  35:  .0625z  <  u  ^  <   .0714z  (A. 14,15) 

w^^^  +  b^^  <   .5  (A. 16) 

35    35 
w^   +  b    <   .5  (A. 17) 

35     35       35 
w^   -  w^   +  150u    =   .5m(3,4)        (A. 18) 

w  "^^  +  b  ^  <   .5  (A.  19) 

35     35       35 
w^   -  w^   +  250u    =   .5m(4,5)        (A. 20) 

artery  56:  .0625z  <  u    <   .07142  (A. 21, 22) 

56    56 
w    +  b    <   .5  (A. 23) 

w,^^  +  b^^  <   .5  (A. 24) 

o         — 

56     56       56 
w    -  w^   +  SOOu    =   .5m(5,6)         (A. 25) 


47 
artery  47:  .0625z  <  u    <   .0714z  (A. 26, 27) 

47    47 
w^   +  b    <   .5  (A.  28) 

47    47  47 

w^   +  b    <  1  -  r^  (A.  29) 

Red  time  ^l^  ^   ^7^^   =   ^  (A. 30) 


constraints 


.4  <  r^^  <   .6  (A. 31, 32) 

25z  <  r^^^  <  50z  (A. 33, 34) 

47    47      47  47 

w,   -  w^  +  500u  =  .5m(4,7)-.5(.5-r  )  (A. 35) 
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artery  16:  .0625z  <  u    <   .0714z  (A. 36, 37) 

16    16  ,   ,„, 

w    +  b    <   .5  (A. 38) 

w/Sb^^  <  1  -  r/^  (A.  39) 

w^^^-  w^^^+  150u^^=  .5m(l,7)-.5(.5-r^^S  (A. 40) 

w^^Sb^^  <   .5  (A.41) 

"7^^'  ''6^^'^  250u^^=  .5m(7,6)-.5(r^^^-.5)  (A. 42) 


cycle  C(l,3,4,7):       m(l,2)  +  m(2,3)  +  m(3,4)  +  m(4,7)  -  m(l,7)  = 

2c(l, 3,4,7)  (A. 43) 

Cycle  C(4,5,6,7):       m(4,5)  +  m(5,6)  -  m(7,6)  -  m(4,7)  = 

2c(4,5,6,7)  (A. 44) 

m(i,j),  c(i,j,k,^)   integers 
",'^  b'J    >  0 

In  using  the  branch  and  bound  algorithm,  the  constraints  were 
introduced  in  the  following  blocks,  each  terminating  with  an  integer 
variable  constraint:    (1)   A.l-A.lO,   A. 16,  A. 23,  A. 28,  A. 38,  A. 11  . 
(2)   A. 12,  A. 13  .   (3)   A. 14,  A. 15,  A. 17,  A. 18  .   (4)   A. 26,  A. 27, 
A.29-A.35  .   (5)   A. 36,  A. 37,  A. 39,  A. 40,  A. 43  .   (6)   A. 19,  A. 20  . 
(7)   A.  41,  A. 42  .   (8)   A. 21,  A. 22,  A. 24,  A. 25,  A. 44  . 
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,44500 


20750 


.35560 


m(l,2) 


m(2,3) 


36405 


m(3,4) 


m(4,7) 


m(l,7) 


35010 


ni(4,5) 


34897 


m(7,6) 


m(5,6) 


Figure  7.   Tree  for  netv.;ork  example.   Labels  on  the  nodes  are  the  calculated 
upper  bounds  on  the  objective  function.   The  objective  function  is  a  weighted 
sum  of  arterial  bandwidths.   Numbers  inside  the  nodes  are  values  for  the  m. 
shown  at  right.   Cycle  constraints  have  been  used  to  limit  the  possible 
choices  of  m(1^7)  and  m(5,6). 


-34 


Appendix  2 

The  mixed-integer  program  for  the  10  signal  problem  is  given 
below. 


Find  b.  z,  w. ,  t  ,  m   to 
1   i   i 


max  b 


Subject  to: 

.0182  <  2  <   .0133  (A2.1,2). 

w^  +  b  <  .53  (A2.3) 

(A2.4,5) 
w^  +  b  <   .60  (A2.6) 

w^  -  w^  +  t^   =   (l/2)m^  -  .035  (A2.7) 


9.4z  <  t   <  12. 5z 


11. 9z  <  t^  <  15. 9z  (A2.8,9) 

-2.03Z  <   .787t2  "  t^  <  2.03z  (A2.10,ll) 

w^  +  b  <   .60  (A2.12) 

"2  "  ^3  ■*"  "2  "   <l/2)m2  (A2.13) 

18. 7z  <  t^  <  25. Oz  (A2.14,15) 

-2.58Z  <   .636t3  -  t^     <  2.58Z  (A2.16,17) 

w^  +  b  <   .53  (A2.18) 

"^3  ■  "4  "*■  S  =   (1/2)^3  +  .035  (A2.19) 
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11, 9z     <     t,      <     15. 9z 

—        4     — 


-4.05Z     <     1.571t.    -   t^     <     4.05Z 
—  ^  3     _ 


w^  +  b     <      .52 


w^   -   w^  +  t^     =      (l/2)m^  +   .005 


(A2.20,21) 
(A2.22_,23) 
(A2.24) 
(A2.25) 


13. 6z     <     t        <     18. 2z 


■2.58z     <      .873t,    -   t,      <     2.58z 
—  5  ^     _ 


w^   +  b     <      ,58 
6  — 


^  ■  ^  ■*■  s 


(l/2)ra^    -    .03 


(A2.26,27) 
(A2.28,29) 
(A2.30) 
(A2.31) 


11. Iz     <     t,      <     14. 8z 
—       0     — 

-2.95z     <     1.231t,    -   t,     <     2,95z 


w     +  b     <      .60 


"6   ■  "7  -^  ^6 


(l/2)ra^   -    ,01 


(A2.32,33) 
(A2.34,35) 

(A2.36) 

(A2,37) 


6.8z     <     t        <     9.1z 


■2.40z     <     1.624t^    -   t^     <     2.40z 

~  10       — 


w„  +  b     <      .60 

o 


W7    -   Wg  +  t^      =      (l/2)m^ 


(A2.38,39) 
(A2,40,41) 
(A2,42) 
(A2,43) 


11.92  <  t-  <  15.92 
o  — 


-1.482  <   .573t„  -  t^  <  1.48z 
—       o     /  — 


w  +  b  <   .60 


"8  ■  ^9  +  'a 


(l/2)m 


8 
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(A2.44,45) 
(A2. 46,47) 
(A2.48) 
(A2.49) 


7.7z  <  tg  <  10.22 


-2.582  <  1.555t_  -  t„  <  2.58z 


w   +  b  <   .58 


"9  '  "lO  "*■  S   =   (^/2)mg  +  .01 


(A2.50,51) 
(A2.52,53) 
(A2.54) 
(A2.55) 


w..  b  >  0 
X  — 


The  above  data  is  arranged  in  blocks.   The  first  (r-1) 

blocks  represent  an  r-signal  problem  for  signals  S  .--..S  . 

1      r 
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